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A four-state model is proposed for description of the sequence of phase tran- 
sitions in ferroelectric crystals of DMAGaS and DMAA1S type. The ordering 
processes in the subsystems of DMA groups are considered as a main rea- 
son of such transformations. The interaction between groups in their various 
orientational states is taken into account in the dipole-dipole approximation. 
Obtained thermodynamical characteristics of the model (spontaneous polariza- 
tion, occupancy of orientational states, dielectric susceptibility, phase diagram) 
are in good agreement with experimental data. The experimental fact of sup- 
pression of ferroelectric phase at increase of hydrostatic pressure is explained 
under assumption that pressure changes the difference between energies of var- 
ious orientational states of DMA groups mostly. 
PACS 77.84.-s, 64.60.Cn 
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1 Introduction 

A family of isomorphous crystals with ferroelectric properties to which e.g. 
(CH 3 ) 2 NH 2 A1(S0 4 )2 • 6H 2 (DMAA1S) and (CH 3 ) 2 NH 2 Ga(S0 4 ) 2 • 6H 2 (DMA- 
GaS) belong to is experimentally studied during the past decade. A peculiar feature 
of the family is a possible existence of the crystal in three different phases at change 
of temperature: at room temperature the crystal is paraelectric, at lowering of 
temperature it sequentially becomes ferroelectric and antiferroelectric (Pietraszko, 
Lukaszewicz and Kirpichnikova, 1993; Pietraszko and Lukaszewicz, 1994; Pietraszko, 
Lukaszewicz and Kirpichnikova, 1995; Tchukvinskyi, Cach and Czapla, 1998). For 
example, it has been measured T& = 136 K and T c2 = 113 K by Tchukvinskyi et al. 
(1998) (or 122 K and 114 K correspondingly (Pietraszko et al, 1995)) for DMAGaS 
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Figure 1: Projection of N and C atoms in the dimethylammonium groups of DMA- 
GaS crystal at ambient temperature (Pietraszko et ai, 1993) onto the ac plane (the 
b coordinate is indicated for each atom; k is the index of DMA sublattices). 



crystal, but it has been only found T c % = 150 K (Pietraszko et ai, 1993; Pietraszko 
et ai, 1994; Pietraszko et ai, 1995) for DMAA1S. The phase transition between the 
ferroelectric and antiferroelectric phases is of the first order (Pietraszko et ai, 1995; 
Tchukvinskyi et ai, 1998). There are evidences, supplied by optical, ultrasonic, 
pyroelectric, dilatometric and dielectric measurements, that the phase transition 
paraelectric - ferroelectric is of the first order close to the second one in DMAGaS 
and of the second order in DMAA1S. Crystallographic analysis shows that in all 
three phases the crystal belongs to monoclinic space groups: the high-temperature 
paraelectric phase has P2i/n space group (Pietraszko et ai, 1993; Pietraszko et 
ai, 1994), ferroelectric and antiferroelectric ones have Pn (Pietraszko et ai, 1994; 
Pietraszko et ai, 1995) and P2i (Pietraszko et ai, 1995) groups respectively. It 
should be mentioned that low-symmetry space groups are subgroups of the high- 
symmetry group obtained by loss of rotation axis and mirror plane respectively 
(point symmetry group 2/m changes to m or 2). 

Certain ferroelectric compounds have structural elements which reorientation 
leads to polarisation of the crystal. Changes in the rest of the crystal are insignificant 
and can be neglected in description of the phenomenon. Crystals DMAA1S and 
DMAGaS belong to this type of ferroelectrics. Here the element, which can be 
reoriented, is the dimethylammonium cation (more strictly NH2 group) (Pietraszko 
et ai, 1995; Sobiestinskas, Grigas, Andreev and Varikash, 1992; Kazimirov et ai, 
1998; Dolinsek et ai, 1999). This element can occupy four equilibrium positions 
related in pairs by inversion centre forming a slightly deformated cross (Figures |l 
and 0). In the paraelectric phase a site in one pair ((k, 1) and (k,2); k = 1,2) is 
occupied with probability 40% and in another ((fc, 3) and 4)) with probability 
10% at 300 K (Pietraszko et ai, 1993). Instead of the presented in Figure p] oblique- 
angled crystallographic coordinates we shall use rectangular ones hereinafter. Their 
OY axes coincides with the b axis, the OX axis is directed along the ferroelectric 
axis (Nl-Nl) and the OZ axis is chosen to make planes XZ and ac identical. 

According to the experimental data (Pietraszko et ai, 1995) in the ferroelec- 
tric and antiferroelectric phases only positions in the first pair are occupied with 
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Figure 2: Projection of N atoms onto the plane perpendicular to the C-C axis of 
the dimethylammonium group for the same crystal as in Figure |l|. 

Table 1: Conversion table of occupancies of orientational states into symmetrized 
ones which correspond to irreducible representations of the point symmetry group 
2/m. 
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appropriate ordering. 

In the ferroelectric phase positions (1,2) and (2,2) (or (1,1) and (2,1)) are mainly 
occupied while in the antiferroelectric state occupancies of positions (1,2) and (2,1) 
(or (1,1) and (2,2) respectively) prevail (Pietraszko et al, 1993). In (Stasyuk and Ve- 
lychko, 2000a) we have performed the symmetry analysis of functions being linear 
combinations of position occupancies of NH2 groups in an elementary cell. Sym- 
metrized linear combinations transforming according to irreducible representations 
of point symmetry group 2/m of the high-temperature phase are presented in Ta- 
ble [TJ. Combinations y + , w__ (?/_, u + ) belonging to representation B u (A u ) correspond 
to orderings appearing in the ferroelectric (antiferroelectric) phase and form order 
parameters for these phases. Existing orderings are of a mixed nature: ferroelec- 
tric ordering in positions (k, 1) and (k, 2) (along the OX axis) is accompanied by 
antiferroelectric one in (k,3) and {k,4) (along the OY axis) and vice versa. 

Proceeding from these facts we have proposed microscopic four-state model (Sta- 
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syuk et al, 2000a) for description of thermodynamics of order- disorder phase tran- 
sitions in the subsystem of reorientationable DMA groups. The phase transition 
between the paraelectric and ferroelectric phases in the considered crystals was de- 
scribed and criteria of order of this transition were established in the approach of 
a direct interaction between variously oriented groups (a site-site approach). The 
results obtained can not be extended on a wide temperature range including e.g. 
the region of the low-temperature phase transition in DMAGaS because total oc- 
cupancies of "longitudinal" positions 1 and 2 (and "transverse" positions 3 and 4 
respectively) have been assumed as temperature independent parameters of theory. 
Nevertheless in this work order parameters corresponding to anti- and ferroelec- 
tric orderings have been constructed what allowed to describe the phase transition 
between ferro- and paraelectric phases and to establish criteria of its order. 

Here a continuation of this work is presented where an approach of dipole-dipole 
interaction between ionic groups is used; temperature change of total occupancies is 
taken into account as well as the effect of reorientational hopping of groups between 
positions (k,p). Ordering of ionic groups corresponding to order parameters of both 
B u and A u symmetry is described in the mean field approximation. Equilibrium 
states are investigated, temperature dependences of order parameters (spontaneous 
polarization) and occupancies of positions of NH 2 groups are studied. Both phase 
transitions (at temperatures T c \ and T c2 ) are described; the regions of existence 
of para-, ferro- and antiferroelectric phases depending on temperature and model 
parameters are established. Observed by experiment changes of picture of phase 
transitions in DMAGaS crystal under the influence of hydrostatic pressure are in- 
terpreted in the framework of microscopic description. 



2 Hamiltonian of the model 



Let us formulate Hamiltonian of the subsystem of DMA groups on the basis of four 
orientational states with taking into account of differences of energies of nonequiv- 
alent configurations, possibility of orientational hopping and interaction between 
groups in dipole approximation. 
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is the Hubbard operator describing transition from the state s' on the state s of 
the DMA complex in the lattice site n and the sublattice k (X*% = \nk, s)(nk, s'|), 
D% k is the a-component of the dipole moment of complex, E a is the corresponding 
component of external electric field, $^'(nn') is the energy of the dipole interaction, 
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E\ and e 2 are the energies of DMA groups in positions (k, 1), (k, 2) and (k, 3), (k, 4) 
respectively, Q\ and Q 2 are the parameters of the reorient at ional hopping. 

Here the coordinate system with the X axis along the axis of spontaneous polar- 
ization of the crystal and the Y axis perpendicular to the DMA plane (parallel to 
the crystallographic b axis) is used. 

It should be mentioned that the Hamiltonian ([!]) simplifies the real structure 
of DMA complexes in DMAGaS crystals, which are oriented not exactly along the 
chosen axes (as assumed in definitions of D" k (§)). Next, in our case the dipole- 
dipole interaction is only the main term of the multipole expansion. Nevertheless, 
as will be shown below, the model described by the Hamiltonian (|l|) can give a 
satisfactory explanation of the behaviour of DMAGaS and DMAAlS-type crystals 
despite of simplifications made. 

The Hamiltonian ([!]) can be rewritten to separate the mean field part 

H = NH C + H MF + H', (4) 
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N is the total number of sites. Further considerations will be restricted to the mean 
field approximation (MFA) and hence the fluctuative term H' is neglected. 

It is convenient to make a unitary transformation to linear combinations of dipole 
moments 

(D n ) = Y,U vM (D a k ) (7) 

ka 

transforming according to the irreducible representations of the high-temperature 
symmetry group of the crystal (cf. Appendix A). After that the MFA Hamiltonian 
looks like 
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E n = U V) kaE a , Ei — E 3 — E 5 — 0, 

feas 

£ 2 = V2E y , £ 4 = v 7 ^*. , E 6 = V2E Z . (11) 

Symmetrized averages of dipole moments form two pairs belonging to different 
irreducible representations: 

A u : (Dx) = ((Df) - (L>I»/V2, (D 2 ) = {(D\) + (Dl))/V2; (12) 
B u : (D 4 > = {(Df) + (£>!»/^2, (As) = {(D\) - (I%))/V2. (13) 

Matrix ip m i has a block-diagonal structure with nonzero blocks belonging to the 
given representation (see Appendix A). 

As follows from the expressions (|1^) and ( |I3"D the antiferroelectric ordering along 
the X axis is accompanied by the ferroelectric one along the Y axis and vice versa. 

In general case / the diagonalization procedure should be applied to the 
matrix in the expression fllPl): 

r 

where 
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(Df) = £ dZ(Xr), ~d r L = £ KtC'Kt. (16) 

r ss' 

Now all necessary thermodynamical functions could be derived from the Hamiltonian 
(H) with corresponding terms expressed by (H) and (|i~4]) . 



and 



3 Thermodynamics of the model 

The partition function of the model in general form is presented as 

Z = Z?Z? exp(-(3NH c ), (3 = l/k B T, (17) 



where 



Z k = £exp(-/3A£), fc = l,2; (18) 

s 

H c = ^(a 1 (D 1 ) 2 + d 1 (D 2 ) 2 + 2b 1 (D 1 }(D 2 ) 

+ a 2 (D 4 } 2 + d 2 (D 5 } 2 + 2b 2 (D 4 )(D 5 }) . (19) 



As a result free energy per particle looks like 



F-~^ = Hc- 1 fnZ 1 - l -^Z 2 . (20) 
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Now one can calculate averages 

(X r k r ) = exp(-/3A£) /j>xp(-/^) 



(21) 



and after substitution into formulae fll2 ), (|T3| ) and (|T6| ) selfconsistency equations for 
(D n ) can be obtained. 

In the case Qi = Q 2 = the matrix Bk is diagonal and expressions for thermody- 
namical functions are more simple. For example, the set of selfconsistency equations 
looks like 
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Components of static dielectric susceptibility are determined in usual way 
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where P x and P y are the components of spontaneous polarization and v c is the 
volume of an elementary cell. After straightforward but cumbersome calculations 
one can derive them from set of equations (^) in the explicit form 
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4 Numerical results 



The model under consideration demonstrates a complicated thermodynamical be- 
haviour and a some numbers of qualitatively different phase sequences dependently 
on parameter values. Further study will be devoted to the case of the sequence of 
phase transitions which is characteristic to the DMAGaS crystal (antiferroelectric 
— > ferroelectric — > paraelectric phases with temperature increase). 

Such a sequence of phase transition is achieved at the next relationships between 
parameters 

< bi, di, a% < b 2l d 2 ; 

ax > a 2 > 0. (26) 

If the above nonequalities are fulfilled, antiparallel orientation of D x - and D^-dipoles 
of neighbour DMA groups is energetically favourable. At low temperatures (when 
occupancies of positions (k, 3) and (k, 4) tend to zero) the antiferroelectric phase 
related with antiparallel dipoles D x is realized. In the result of increased probabil- 
ity of appearance of groups in the mentioned positions at temperature growing the 
interaction with participation of D y dipoles becomes significant. The phase tran- 
sition to the phase with their antiferroelectric orientation occurs when respective 
interactions become prevailing. This phase is simultaneously ferroelectric along the 
OX axis. 

There are two other more subtle conditions: for a correct description of the 
DMAGaS crystal the high-temperature phase transition should be of the first order 
and the ratio (6 c i — C 2)/@ci should correspond to the experimental one. 

Results of numerical calculations of temperature dependences of symmetrized 



dipole moments ([12]) , ( |l3l) (obtained as solutions of equation set (|3)) having mean- 
ing of order parameters are depicted in Figure |3|. Here the case of the first order 
high-temperature phase transition close to the second one is illustrated. Besides 
the temperatures of real phase transitions c i and Q c2 , the temperature of a possi- 
ble phase transition 0^, (which might be realized in the absence of the ferroelectric 
phase) is indicated. 

As one can see in Figure |3] and from the expressions flT2|), (|I~3|) antiferroelectric 
orientation along the X axis is accompanied by ferroelectric one along the Y axis and 
vice versa. Polarization along the X axis is larger than along the Y one due to larger 
occupancies of positions 1 and 2 comparatively to positions 3 and 4 (Figure §). 

Behaviour of the spontaneous polarisation P x (which is proportional to (-D4) 
as stated in (p3|)) strictly corresponds to experimental measurements (Dacko and 
Czapla, 1996; Pykacz and Czapla, 1997). But despite of common assumption (see 
e.g. Kazimirov et al, 1998) at lowering of temperature the occupancy of the positions 
3 and 4 do not tend to zero monotonically: occupancies are asymmetrical and (Xf A ) 
even increases in the ferroelectric phase (with vice versa occupancies for k = 2). 

The influence of the orientational hopping on the phase transition in the system 
resemble one for Ising model (Figure |5] - for simplicity is assumed Qi — Q 2 — 
Namely the increase of Q value leads as a whole to the decrease of temperatures 
of both phase transitions and their sequential zeroing at corresponding values fl c2 
and Q c \. Antiferroelectric phase vanishes at relatively low values of Q when the 
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Figure 3: Dependences of (D) on temperature at specified below values of param- 
eters: &i/ai=1.55, di/ai=0.5, a2/ai=0.4, b 2 /ai=2.6, d 2 /ai=1.6, A/aic^=i.5855, 
r2i=Sli=0, d y /d x = 1.4; 9 is given in a\d\ units. The same values of parameters 
are held in Figures Solid lines correspond to the thermodynamically stable 

solution, dashed and dotted - to the unstable ones. 
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Figure 4: Dependences of occupancies of states on temperature. Here k = 1. 
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Figure 5: Phase diagram 6 — VI. and Figure 6: Phase diagram 6 — A. and 
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temperature of the other phase transition changes insignificantly. Due to prevail- 
ing suppression of antiferroelectric ordering and the mentioned above competition 
between anti- and ferroelectric interactions the considered diagram has a peculiar 
feature: the temperature of the ferroelectric-paraelectric phase transition increases 
at small value of Q. 

According to the experimental data on pressure effects in DMAGaS (Yasuda, 
Kaneda and Czapla, 1999) applying and increase of hydrostatic pressure leads to 
the narrowing of the temperature region of existence of the ferroelectric phase. At 
certain pressure the ferroelectric phase is completely suppressed. Thermodynamical 
description in the framework of the Landau theory allows to reproduce observed 
phase diagram by setting coefficients at the second order terms in the Landau ex- 
pansion of free energy to be pressure dependent (Stasyuk, Velychko, Czapla and 
Czukwinski, 2000b). Microscopic description of influence of hydrostatic pressure is 
rather a complicated question. One can assume, that due to limited compressibility 
of the crystal interaction constants change insignificantly and the main effect is due 
to the increase of the difference of site energies A. In this case the obtained phase 
diagram remarkably good coincides with qualitative form of the experimental one 
(Figure ^). All three phases could exist at change of temperature in the region 
bounded by values A cl and A c2 . This confirms the above assumption of prevailing 
influence of pressure on site energies. Nevertheless, it seems that a more complete 
reproduction of the observed dependences (in particular, the increase of temperature 
of phase transition in the antiferroelectric phase at the pressure values above the 
triple point) can be achieved when additional factors (such as change of interaction 
constants or the parameter Q) will be taken into account. For example, values of 
parameters of interaction are known to increase at rising of pressure. The tempera- 
ture presented in Figures |3|-|8] is made dimensionless by division on a\d 2 x . Hence if 
the dimensionless C remains constant at increase of pressure, the real temperature 
of the phase transition grows up. 

As one can mention there are both triple and tricritical points on the diagram. 
There is a narrow region of A values near the triple point where both phase transition 
are of the first order as in the case of DMAGaS crystals. 

Parameters in Figures |3] and |] are chosen to fulfill the mentioned above conditions 
and to give the best description of experimental data for DMAGaS crystal. On the 
other hand, a typical behaviour of the DMAA1S crystal could be obtained in two 
different way. One can assume that value of the parameter A is smaller than A C 2 
for this case (Figure |6[). Alternatively, value of the parameter Q can be less than 
Q C 2 (Figure |5|). In both cases it gives direct transition of the second order from the 
paraelectric phase to the ferroelectric one at lowering of temperature. 

Dependences of inverse components of dielectric susceptibility are presented in 
Figure [T[ Xx(T) has a pronounced peak at the temperature of phase transition 
between ferroelectric and paraelectric phases but does not reach singularity. Such 
a behaviour good coincides with the experimental measurements of dielectric per- 
mittivity made by Tchukvinsky et al. (1998). Due to the first order of the high- 
temperature phase transition the slope of the \ x (T) curve changes more than twice 
in the phase transition point. Xy{T) has a small peak in this point (Figure [8]). There 
is a lack of experimantal data considering dielectric measurements along the OY axis 
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Figure 8: Dependence of the Y component of static dielectrical susceptibility on 
temperature. 



for the DMAGaS crystal but those for the DMAA1S (Kapustianik et al, 1994) close 
resemble the theoretical picture. Both components have jumps at temperature of 
the first order phase transition from antiferroelectric to ferroelectric phase. 

Basing on the experimental data some approximate evaluation of values of pa- 
rameters has been made. A value of polarization of saturation Pg = 1.8 ■ 10 -2 C/m 2 
gives d x = 6.8 • 10~ 30 C-m (what corresponds to 0.9e x 0.5 A). Taking an approx- 
imate value ai = 27 'A ■ 10 20 cm -3 , we obtain aid 2 x = 8.3 ■ 10~ 3 eV (normalizing 
factor for 6, A and ft) and a x v c = 27 'A ■ 10 20 cm~ 3 x 786 • 1(T 24 cm 3 = 2.15 (nor- 
malizing factor for the components of x)- Hence we get A = 0.013 eV (this value 
corresponds to the observed ratio of occupancies of positions 1,2 and 3,4 at ambient 
temperature), VL c2 — 40 cm -1 and the maximum value of Xy = 13 at the point of 
the high-temperature phase transition what good coincides with the experiment for 
DMAA1S. 
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5 Conclusions 



A simple dipole approximation for the four-state model proposed for description 
of thermodynamics of the DMAGaS-DMAAlS family of crystals is developed. The 
model exhibits complicated thermodynamical behaviour depending on values of sys- 
tem parameters. At certain values of parameters the model gives the same sequences 
of phase transitions as found experimentally in DMAGaS and DMAA1S. 

Numerical exploration of the model shows that increase of the orientational hop- 
ping between sites leads to the suppression of ordered phases in the Ising model like 
way. 

An attempt to give theoretical explanation of the experiment on the influence of 
hydrostatic pressure on phase diagram of DMAGaS crystal results in good qualita- 
tive agreement. It is assumed that the main effect of pressure is the change of the 
difference of site energies. 
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7 Appendix A. Transformation of the matrix of 
dipole-dipole interactions to irreducible repre- 
sentations of the high-temperature symmetry 
group 
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Here axis indices are changed first (e.g. {ka} = {lx, ly, lz, 2x, 2y, 2z}). 
Matrix of unitary transformation is as follows: 
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After the transformation the matrix of interactions 



1> 



( 


a'-d' 


b' + e' 


c' 






b' + e' 


g' + k' 


/i' 


-r 




c'-f 


h'-r 


mf 


— n' 




























\ 













( 


di bi 


ci 





\ 




bi di 


ei 










ci ei 


/i o 













a 2 


62 


c 2 







6 2 


d 2 


e 2 


\ 





c 2 


e 2 


h) 



\ 





a' + d' b'-e' c' + f 

b'-e' g'-k' h' + V 

c' + f h' + 1' n' + m' J 



(29) 



becomes block-diagonal, where the upper left block corresponds to the A u represen- 
tation and the bottom right one to the B u . 
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